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ABSTRACT 

This  thesis  develops  a  technique  for  the  three-dimensional  (3-D)  display  of 
clouds,  topography  and  atmospheric  surfaces  on  an  image  and  graphics  analysis 
system.  Cloud  heights  were  estimated  from  infrared  satellite  imagery.  Software  was 
written  to  translate  two-dimensional  gridded  data  height  fields  into  gray  shade  images. 
Existing  3-D  display  software  was  modified  to  use  the  data  images  as  input.  Potential 
temperature  and  500  mb  surfaces  are  presented  as  examples  of  possible  fields  for  anal- 
ysis in  3-D.  The  software  was  further  modified  to  expand  the  geographical  area 
enclosed  in  the  final  3-D  images.  These  images  were  used  to  produce  3-D  displays  of 
topography,  clouds  and  atmospheric  parameters.  The  relationships  between  atmos- 
pheric surfaces  and  cloud  structures  can  be  more  clearly  seen  in  the  3-D  representa- 
tions. With  modifications  to  increase  the  resolution  of  the  3-D  displays,  improved 
understanding  of  the  atmosphere  can  be  expected. 
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I.  INTRODUCTION 

During  the  past  40  years,  significant  advances  have  been  made  in  the  modeling 
and  understanding  of  atmospheric  processes.  Numerical  techniques  have  advanced  at  a 
pace  closely  linked  to  the  development  of  more  sophisticated  computer  systems 
(Shuman,  1978).  Remote  sensing  of  the  atmosphere  has  the  potential  to  further  the 
accuracy  of  numerical  models  and  meteorological  prognoses  by  providing  more 
complete  data  bases.  However,  one  area  that  does  not  appear  to  have  advanced  in  like 
manner  is  the  graphical  and  visual  modeling  of  the  atmosphere. 

For  most  of  the  past  four  decades,  meteorologists  have  depended  on  two- 
dimensional  representations  of  atmospheric  variables.  Surface  and  atmospheric  data 
are  plotted  either  by  computer  or  manually  on  varying  sizes  of  surface  maps.  On  each 
map  is  overlayed  wind,  pressure  and/or  temperature  information  at  some  vertical  level. 
This  may  be  of  observe.  ^recast  (predicted)  patterns.   The  individual  meteorologist 

must  then  keep  track  of  these  various  layers  by  mentally  overlaying  them.  If  the  size 
of  the  maps  agree,  it  is  sometimes  possible  to  physically  overlay  transparent  versions. 
This  stacking  process  leads  the  educated  eye  to  make  associations  between  the  changes 
at  various  levels. 

A  surface  front,  for  example,  may  be  strongly  associated  with  an  upper-level 
trough.  A  developing  area  of  severe  thunderstorms  is  usually  the  result  of  a  very 
complex  interaction  between  dry  and  moist  air  low  in  the  atmosphere,  the  jet  stream, 
surface  heating  and  upper-level  troughs.  The  numerical  models  which  generate  these 
charts  have  improved,  but  their  output  is  still  in  the  two-dimensional  (2-D)  map  form. 

The  need  for  these  models  is  clear  when  instructing  students.  The  models  are 
used  to  assist  students  in  conceptualizing  the  3-D  nature  of  the  atmosphere  and  the 
interrelationships  of  its  constituents.  Three-dimensional  (3-D)  constructions  have  been 
limited  mostly  to  tedious,  artistic  renderings  produced  by  hand,  or  physical  models 
constructed  from  plastics.  Advanced  3-D  computer  techniques  for  representing  the 
environment  have  become  available  only  in  the  past  10  years. 

Attempts  have  been  made  to  use  stereo  techniques  to  assist  users  in  perceiving  a 
more  complete  representation  of  the  environment  (Hasler  et  al,  1981).  In  the  stereo 
technique  the  viewer  is  looking  at  two  images  simultaneously:  the  original  image  in  red 
and  a  shifted  computer  generated  blue-green  version  of  the  image.    In  the  shifted 


image,  each  picture  element  (pixel)  pixel  is  moved  to  the  right  of  its  real  position  by  a 
distance  that  will  correspond  to  the  height  changes  (depth)  in  the  displayed  field.  The 
viewer  wears  colored  lenses  to  properly  sense  the  depth  in  the  field.  A  blue-green  lens 
covers  the  left  eye  and  a  red  lens  covers  the  right.  Each  eye  then  sees  one  of  the  two 
images,  the  "real"  pixel  positions,  and  a  shifted  pixel  position.  The  result  is  a  percep- 
tion of  depth. 

Three  serious  limiting  factors  affect  the  use  of  this  type  of  product.  First,  only  a 
very  small  part  of  the  total  satellite  data  base  is  processed  for  stereoscopic  display. 
The  user  is  not  free  to  construct  and  view  pictures  of  his  or  her  designation.  Second, 
the  user  must  wear  special  glasses  to  view  the  images.  Even  with  these  lenses,  it  still 
requires  some  interpretive  skill  on  the  part  of  the  user  to  achieve  the  proper  sense  of 
depth  in  all  representations.  Third,  the  product  is  static  and  cannot  be  manipulated  by 
the  viewer.  Inherent  in  research  and  operational  activities  is  the  requirement  that  the 
product  be  as  dynamic  as  possible;  that  it  can  be  manipulated  with  as  much  speed  as 
possible  to  test  various  hypotheses. 

The  products  used  in  meteorology  require  the  user  to  establish  the  proper  3-D 
relationships  from  2-D  products.  By  using  techniques  originally  designed  for  mechan- 
ical engineers  and  architects,  computer  displays  can  be  produced  which  place  atmos- 
pheric data  in  a  3-D  perspective. 

A  first  step  in  this  direction  for  meteorology  was  taken  at  Colorado  State 
University  (CSU)  (Meade,  1985).  Based  on  a  3-D  computer  display  program 
MOVIE. BYU,  originally  developed  at  Brigham  Young  University,  subroutines  were 
added  which  allowed  the  construction  of  models  of  clouds  over  a  topographic  surface. 
After  modification,  it  was  renamed  MOVIE2.  The  input  was  a  combination  of  infrared 
and  visual  satellite  data  from  GOES  satellites,  plus  a  digitized  topography  of  the  conti- 
nental U.S.,  from  the  Defense  Mapping  Agency  (DMA).  The  infrared  imagery 
provided  cloud  top  height  information  while  the  visual  data  were  used  optionally  to 
shade  the  resulting  3-D  structures-.  The  program,  as  it  existed  then,  allowed  only  the 
display  of  a  4.7  by  4.7  degree  latitude  area. 

This  areal  limitation  has  been  removed  here.  3-D  displays  can  now  be  produced 
for  areas  from  as  small  to  as  large  as  desired.  However,  if  the  full  resolution  of  the 
data  for  clouds,  topography  and  atmospheric  data  surfaces  were  used,  the  run  times  for 
production  of  the  3-D  images  would  be  significant.  The  method  used  to  overcome  this 
problem  will  involve  a  simple  averaging  of  picture  elements.   This  process  will  increase 


the  number  of  pixels  averaged  together  as  the  size  of  the  selected  area  increases.  This 
will  reduce  the  effective  resolution  so  that  it  holds  the  image  content  to  a  size  which 
the  MOVIE2  software  will  manipulate  in  a  reasonable  amount  of  time. 

Utilizing  the  capability  of  the  GEMPAK  software  ,  developed  at  the  Goddard 
Space  Flight  Center  and  residing  on  the  NPS  VAX  1 1/780  computer,  standard  vertical 
weather  measurements  (rawinsonde  data)  can  be  plotted  and  gridded.  These  gridded 
data  sets  can  be  of  either  observed  or  derived  quantities  such  as  the  lifted  condensation 
level,  stability  indices,  constant  pressure  surfaces  and  constant  potential  temperature 
surfaces.  Once  gridded,  these  data  should  then  be  transformed  into  a  continuous  tone 
gray-shade  image.  The  image  then  can  be  displayed  on  graphics  devices  similar  to  the 
COMTAL  Vision  1  at  NPS.  This  image  can  be  manipulated  and  displayed,  simultane- 
ously with  cloud  structures,  by  using  MOVIE2  software.  The  combined  display  is 
expected  to  greatly  assist  researchers  in  viewing  large  mesoscale  and  small  synoptic- 
scale  systems.  Further,  it  should  enhance  understanding  of  the  relationships  between 
various  atmospheric  surfaces. 

It  is  the  intent  of  this  thesis  to  make  a  second  step  towards  the  development  of  a 
tool  which  research,  and  eventually  the  general  meteorological  community,  can  use. 
The  technique  developed  in  this  research  will  allow  for  the  near  real  time  display  of 
integrated  topography,  clouds  and  other  atmospheric  information  in  3-D.  These, 
combined  with  conventional  data  display  techniques,  will  lead  users  to  quickly  extract 
information  about  the  state  of  the  atmosphere. 
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II.  SOFTWARE 

There  are  four  main  groups  of  software  used  in  the  construction  of  the  3-D 
models:  GEMPAK,  NAVGEN,  GRIDIMG,  and  MOVIE2.  These  programs  and 
program  packages  take  raw  satellite,  rawinsonde  and  mercator  projected  topography 
data,  and  combine  them  into  a  unified  display  of  clouds,  topography  and  atmospheric 
variable  surfaces.  Although  the  manner  in  which  the  final  3-D  image  is  produced  is 
cumbersome,  future  refinements  to  this  sequence  are  expected  to  greatly  simplify  the 
procedure.  The  flow  diagram  (Fig.  2.1)  traces  the  processing  of  data  through  input  to 
the  3-D  display  software. 

A.      GEMPAK 

GEMPAK  is  a  software  package  developed  at  the  NASA/Goddard  Space  Flight 
Center  (desJardins  et  ai,  1981).  It  consists  of  a  variety  of  utility  programs  designed  to 
operate  on  either  surface  or  vertical  sounding  (rawinsonde)  data  files.  It  allows  the 
user  to  create  two-dimensional  graphical  displays  of  reported  or  derived  atmospheric 
variables.  It  is  menu  driven,  provides  help  files,  and  has  a  complete  instruction  manual 
for  reference.  Only  a  small  part  of  its  capabilities  were  used  in  preparing  data  for 
eventual  3-D  display.  Four  GEMPAK  programs  were  used:  OAGRID,  SNBNANL, 
GDLIST,  and  REMAP. 

1.  OAGRID 

The  program  OAGRID  produces  a  gridded  data  file  to  accept  output  later 
from  the  program  SNBNANL.  OAGRID  reads  a  specified  input  vertical  sounding  or 
surfaces  data  set,  or  a  subset  specified  by  the  user.  It  produces  an  evenly  spaced, 
gridded  field  for  one  or  more  parameters  (e.g.,  heights,  temperatures,  lifted  condensa- 
tion levels,  etc.). 

2.  SNBNANL 

Having  established  a  suitable  output  file  with  OAGRID,  the  actual  objective 
analysis  is  done  by  the  program  SNBNANL.  It  performs  a  Barnes  objective  analysis 
on  the  input  data  (Barnes,  1964).  This  produces  data  at  evenly  spaced  locations.  This 
step  is  necessary  since  rawinsonde  balloon  sounding  launch  sites  are  not  located  at  grid 
points. 
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3.  GDLIST 

This  program  allows  the  user  to  copy  the  data  output  from  SNBNANL  to 
another  disk  file,  the  terminal  screen,  or  to  a  printer.  The  user  can  have  any  or  all  of 
the  gridded  parameters  created  in  SNBNANL  sent  to  these  devices. 

4.  REMAP 

The  GEMPAK  program  REMAP  converts  512  X  512  pixel  satellite  or  atmos- 
pheric data  images  into  various  types  of  map  projections.  The  image,  representing 
infrared  radiance  values,  is  input  with  the  corresponding  navigation  files;  one  file  each 
for  the  latitude  and  longitude  of  image  pixels.  It  is  a  menu  driven  program  prompting 
the  user  at  each  step  for  input  and  output  file  names,  projection  type,  and  borders  (if 
any).   The  resulting  image  is  then  availble  for  use  as  input  to  the  MOVIE2  programs. 

B.  GRIDIMG 

The  next  step  was  the  development  of  GRIDIMG  which  uses  the  gridded  data 
file  as  input  and  performs  interpolation  between  adjacent  grid  points  in  both  the  x  and 
y  direction  (bilinear  interpolation).  The  interpolation  allows  a  continuous  tone  gray 
shade  image  to  be  produced.  Each  shade  represents  a  different  height  value  where  the 
atmospheric  variable  occurs.  In  order  to  assign  the  proper  range  of  gray  shade  values 
and  prepare  the  gridded  data  file,  the  user  must  examine  the  maximum  and  minimum 
heights  in  the  data  file.  The  maximum  and  minimum  gray  shade  must  be  entered  to 
provide  a  range  of  values  to  assign  to  grid  point  heights.  Also,  the  grid  dimensions  (17 
X  17,  12  X  12,  etc.)  and  number  of  times  the  interpolation  is  performed,  per  image  line, 
must  be  added.  These  values  are  then  converted  into  the  IR  radiance  count  values 
associated  with  those  heights,  on  a  representative  rawinsonde  sounding. 

C.  NAVGEN 

Next,  NAVGEN  was  written  so  that  output  gray  shade  images  from  GRIDIMG 
would  have  navigation  files.  These  files  allow  conversion,  through  REMAP,  of  the 
image  into  mercator  coordinates.  It  creates  the  navigation  files  based  on  the  input 
dimensions  of  the  area  contained  in  the  output  file  from  GRIDIMG.  It  produces  two 
output  files  describing  the  latitude  and  longitude  of  every  other  pixel  on  a  line  of  the 
image.  This  is  done  by  matching  the  data  point  heights  to  the  temperature  at  that 
height  on  a  rawinsonde  sounding.  This  temperature  is  converted  into  an  IR  radiance 
value  (gray  shade)  using  a  standard  GOES  satellite  calibration  table. 
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D.      MOVIE2 

The  MOVIE2  software  allows  the  user  to  create  3-D,  mathematically  defined 
structures  on  a  computer  (e.g.  VAX  1 1/780)  for  display  on  a  computer  graphics  device 
(e.g.  COMTAL,  Tektronics,  PC,  etc.).  It  is  a  complex  general  application  program 
allowing  a  user  to  create  3-D  images  by  inserting  the  parameters  necessary  to  define 
the  dimensions  of  a  structure.  The  structure's  computer  representation  then  can  be 
manipulated  by  issuing  specific  commands.  The  shape  can  be  rotated,  enlarged, 
colored  and  separted  into  its  various  components. 

The  CSU  developed  user  software  includes  a  variety  of  programs  which  allow  the 
translation  of  continuous  tone  images,  in  256  shades  of  gray  from  black  (shade  0) 
through  white  (shade  256).  The  CSU  user  programs  DMADAT,  TOPGEN,  DEFSCT, 
INTCLD,  CLDGEN  and  MERGER  were  added  to  allow  the  display  of  digital  satel- 
lite data  rather  than  engineering  or  architectural  data,  for  which  it  was  first  written 
(Meade,  1985).  DISPLAY,  UTILITY,  SECTION,  TITLE,  MOSAIC,  COMPOSE  and 
UPDATE  are  the  original  FORTRAN  programs  which  made  up  the  original 
MOVIE. BYU  software  collection.  As  DISPLAY  is  the  main  program  and  can  be  used 
independent  of  the  others,  this  program  was  used  to  display  the  3-D  images. 

1.  DMADAT 

This  program  allows  the  user  to  select  a  subarea  of  the  data  contained  on  the 
DMA  topography  data  tape.  Based  on  the  user's  answers  to  prompts  concerning  the 
location  and  size  of  the  desired  area,  DMADAT  reads  in  the  appropriate  tape  records. 
The  input  data  are  converted  into  gray  shades,  corresponding  to  temperature  and 
vertical  height,  and  writes  the  output  image  data  to  a  scaled  topography  file.  At 
present,  the  area  selected  can  range  from  1  X  1  degree  up  to  20  X  20  degree  areas. 

2.  TOPGEN 

The  DMADAT  output  file  must  be  converted  into  a  form  useable  by  the  3-D 
DISPLAY  software.  TOPGEN  performs  this  function  by  converting  the  input  topog- 
raphy image  into  a  mathematical  description  of  the  data.  It  requires  that  the  user 
enter  the  portion  of  the  topography  image  to  be  included  in  the  DISPLAY  readable 
file.  This  can  be  either  the  entire  image  or  a  subset.  The  number  of  polygons  desired, 
on  a  side  of  the  3-D  structure,  and  a  vertical  height  exaggeration  factor  are  entered. 
The  output  file  may  be  displayed  either  by  itself,  or  merged  with  clouds  and  atmos- 
pheric surfaces. 
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3.  DEFSCT 

The  program  DEFSCT  keeps  track  of  which  files  have  been  requested  for  use 
in  the  final  merged  cloud  and  topography  model.  It  contains  information  on  the  side 
length  selections  for  clouds,  topographies  and  atmospheric  surfaces  as  well  as  the  name 
of  files  to  be  used  for  input  and  output.  Before  execution  of  this  program,  the  output 
file  DEFSCT.DAT  must  be  deleted  so  that  the  program  does  not  maintain  values 
selected  by  the  user  describing  other  previous  3-D  image  constructions  (side  lengths 
and  file  names).  DEFSCT  is  executed  after  TOPGEN  and  before  any  other  operations 
are  performed  on  a  cloud  or  upper  air  data  image  file  (Meade,  1985). 

4.  INTCLD 

Cloud  top  heights  can  be  determined  from  the  infrared  (IR)  temperature  of 
cloud  tops  by  converting  the  gray  shade  to  the  corresponding  temperature.  The 
temperature  is  then  converted  into  a  height  by  comparing  it  to  the  heights  and  associ- 
ated temperatures  on  a  representative  rawinsonde  sounding.  These  z-axis  heights  are 
included  in  the  file  to  be  read. later  by  the  MOVIE2  software.  The  x  and  y-axis  extent 
of  a  cloud  are  determined  by  a  software  comparison  of  the  cloud/no-cloud  threshold 
value  set  by  the  user. 

In  the  program  INTCLD  the  cloud/no  cloud  threshold  sets  the  value  at  which 
the  software  should  begin  to  consider  pixel  gray-shades  to  represent  clouds.  It  is 
entered  in  response  to  queries  by  the  CSU  program  INTCLD  (Meade,  1985).  When  a 
meteorological  variable  surface  is  to  be  produced,  a  value  of -999  is  entered  in  response 
to  the  question  of  base  height.  This  entry  sets  the  base  and  top  of  the  layer  at  the 
same  level,  resulting  in  a  curved  plane  having  no  vertical  thickness.  This  program 
reads  each  line  of  the  input  image  file  and  averages  the  pixel  values  together  to  produce 
smoothed  image  files.  A  limit  of  three  levels  of  clouds  and/or  atmospheric  surfaces 
currently  exists  in  the  software.  Two  files  for  each  cloud  level  are  produced  storing 
information  on  the  base  level  of  the  image  and  horizontal  extent  within  the  cloud  or 
surface. 

5.  CLDGEN 

The  program  CLDGEN  produces  the  files  that  can  be  merged  together  with 
topography  files  from  TOPGEN  and  can  be  read  by  DISPLAY.  In  creating  a  meteo- 
rological upper  air  surface,  rather  than  clouds,  the  DEFSCT.DAT  file  must  first  be 
deleted.  This  is  necessary  because  CLDGEN  will  attempt  to  search  for  files  that  have 
not  been  created  in  INTCLD  when  using  the  -999  cloud  base  option.    The  inputs  for 
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this  program  are  the  output  files  from  INTCLD.    The  user  provides  the  names  for 
these  files,  CLDGEN  reads  the  input  data  files  and  produces  the  output  scaled  cloud 
or  upper-air  surface  files. 
6.   MERGER 

The  program  MERGER  combines  the  topography  and  cloud  or  surface  files 
into  a  single  "cloud  topography  model"  file.  If  an  upper-air  surface  is  merged  with 
topography,  the  DEFSCT.DAT  file  must  be  deleted  again,  as  it  is  recreated  at  each 
processing  step  in  the  CSU  added  programs.  Once  created,  the  output  file  can  be  used 
as  input  to  the  program  DISPLAY. 

When  the  final  model  includes  clouds,  topography  and  an  upper-air  surface, 
DEFSCT  and  INTCLD  are  run  separately  for  each  layer.  The  upper-air  surface 
description  replaces  one  pair  of  files  describing  an  unused  cloud  layer  in  the  cloud 
model  files.   When  iMERGER  is  run  the  result  is  a  combined  representation. 

E.       DISPLAY 

DISPLAY  is  a  general  purpose  computer  graphics  display  program  for 
constructing  polygonal  representation  of  objects.  It  has  the  capability  to  produce 
either  line  or  continuous  tone  image  output  (Christiansen  and  Stephenson,  1981a). 
There  are  47  interactive  commands  which  allow  for  changes  to  be  applied  to  a  previ- 
ously created  image  or  frame. 

1 .    Commands 

A  set  of  3-D  images,  once  displayed  by  the  program  DISPLAY,  is  referred  to 
as  a  frame.  Within  that  frame  each  component  is  a  numbered  part  (Christiansen  and 
Stephenson,  1981a).  By  identifying  an  individual  part  of  the  frame,  it  can  be  manipu- 
lated. For  instance,  a  frame  may  consist  of  three  parts;  the  topography,  a  cloud  and 
an  atmospheric  data  surface.  Any  additional  cloud  or  surface  would  be  considered 
another  part  of  the  frame.  Some  of  the  47  commands  can  be  used  to  change  the  entire 
frame  (colors,  apparent  distance,  viewing  angle)  or  may  apply  changes  only  to  specified 
parts  (Christiansen  and  Stephenson,  1981b).  The  ROTATE  command  rotates  all  the 
parts  within  the  frame  by  a  specified  number  of  degrees  on  a  given  axis.  DISPLAY 
does  this  while  maintaining  the  proper  relationship  between  the  components  of  the 
frame,  just  as  if  the  observer  were  to  change  his  or  her  viewing  position.  When  several 
parts  are  introduced  into  a  frame,  coloring  individual  parts  becomes  a  useful  aid  in 
recognizing  them  after  their  positions  have  changed. 
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2.    Geometry  Files 

DISPLAY  can  use  four  types  of  data  files:  geometry,  displacement,  scalar 
function  and  contour.  For  purposes  of  this  thesis,  only  geometry  files  were  used.  Each 
geometry  file  must  contain  a  complete  description  of  the  objects  to  be  handled  by 
DISPLAY.   The  description  includes: 

•  Control  variables  representing  the  number  of  parts  that  will  appear  in  the  frame, 
the  number  of  connection  points  (nodes)  between  line  segments,  and  the  total 
number  of  edees  in  all  the  individual  parts.  A  square  plane  consisting  of  none 
small,  four-sided  panels  would  have  36  edges 

•  A  connectivity  array  associating  the  nodes  with  each  small  element,  or  panel,  of 
which  they  are  a  part 

•  The  coordinate  array  storing  the  location  of  the  panel  connection  points  (nodes) 
in  Cartesian  coordinates 

•  The  parts  array,  a  two-dimensional  matrix  containing  the  upper  and  lower  limits 
of  the  connected  panels  in  a  part. 
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III.  CASE  DATA  AND  PROCESSING 

Selection  of  the  case  to  demonstrate  synoptic-scale  weather  features  in  3-D,  was 
based  on  a  review  of  weather  satellite  images  during  the  period  15  -  17  May  1986. 
Also,  the  availability  of  the  corresponding  rawinsonde  data  played  a  large  part  in  the 
case  selection.  As  discussed  above,  three  primary  types  of  data  are  required:  satellite, 
rawinsonde  and  topography. 

A.  SATELLITE  DATA 

An  infrared  satellite  image  from  the  GOES  satellite  at  2300  hours  Greenwich 
Mean  Time  (GMT)  on  15  May  1986  was  used.  The  area  selected  includes  the  south- 
central  U.S.  in  an  area  bounded  by  25°  and  45°  N  latitude,  and  by  80°  and  100°  W 
longitude.  The  infrared  satellite  image  on  that  date  indicated  well  developed  cloud 
structures,  with  good  vertical  extent,  in  a  cyclonic  weather  system.  The  image  was 
provided  by  the  Naval  Environmental  Prediction  Research  Facility  (NEPRF), 
Monterey,  California.  The  images  were  of  4  n  mi  resolution  at  the  subsatellite  point. 
This  enlarges  to  approximately  6  n  mi  resolution  at  40°  N.  That  is,  the  smallest 
picture  element  (pixel)  distinguishable  represents  a  6  n  mi  square  of  the  earth.  The 
average  radiance  from  the  pixel  area  is  represented  as  a  single  brightness  value. 

Included  with  the  image  were  two  navigation  files.  The  navigation  files  contained 
the  information  necessary  for  a  computer  program  to  associate  each  pixel  with  the 
latitude  and  longitude  of  its  location.  The  raw  image  was  in  the  satellite's  coordinate 
system.  Therefore,  this  could  be  converted  into  a  mercator  projection  for  use  in  the 
MOVIE2  manipulations. 

B.  RAWINSONDE  DATA 

Computer  files  containing  vertical  atmospheric  sounding  data  for  0000  GMT  16 
May  1986,  over  the  continental  U.S.,  were  used  as  inputs  to  the  programs  GDFILE, 
SNBNANL,  and  GDLIST.  These  programs  produced  evenly  spaced,  17  X  17  point 
(60  n  mi  spacing)  gridded  data  files  containing  500  millibar  (mb)  pressure  and  constant 
potential  temperature  (330°  K  and  320°  K)  surfaces.  The  500  mb  pressure  surface, 
occuring  at  roughly  18,000  feet,  is  a  standard  level  used  by  meteorologists  to  determine 
the  location  and  intensity  of  significant  weather  features.  The  constant  potential 
temperature  surfaces  can  also  provide  information  on  the  stability  of  the  atmosphere. 
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The  gridded  data  files  contain  tables  of  numbers  giving  the  height  at  which  the 
given  surface  occurs  in  the  atmosphere.  MOVIE2  requires  gray  shade  continuous  tone 
images  as  input.  To  create  the  continuous  tone  image  from  these  data  files  they  were 
processed  by  the  program  GRIDIMG.  The  output  of  GRIDIMG  is  a  512  by  512 
pixel  continuous  tone  image. 

The  data  images  output  by  GRIDIMG  were  in  a  cylindrical  map  projection.  The 
program  NAVGEN  was  used  to  create  the  navigation  file's  for  the  images.  The  navi- 
gation files  and  image  were  used  as  input  to  the  GEMPAK  program  REMAP  to 
produce  mercator  projected  images. 

C.  TOPOGRAPHY  DATA 

The  topography  data  files  were  obtained  from  the  Defense  Mapping  Agency 
(DMA).  These  files  provided  elevation  data  over  the  continental  U.S.  in  30  second 
(1/120  degree)  intervals.  These  data  were  read  by  the  MOVIE2  program  DMADAT 
(Meade,  1985).  A  subset  of  the  data  selected  was  chosen  to  coincide  with  the 
remapped  satellite  and  atmospheric  data  images.  The  output  from  DMADAT  was  a 
gray  shade  image  in  a  mercator  projection.  This  image  was  used  as  input  to  the 
program  TOPGEN,  which  translated  the  image  into  a  3-D  mathematical  description. 
The  output  from  TOPGEN  was  available  for  immediate  use  by  DISPLAY,  or  for 
merging  with  cloud  and  parameter  surfaces. 

D.  DATA  PROCESSING 

1.  Satellite  Imagery 

The  raw,  infrared  satellite  image  (Fig.  3.1)  is  shown  on  the  COMTAL  image 
processor.  This  initial  image  enclosed  an  area  much  larger  than  was  desired.  The  raw 
image  was  cropped  down  to  a  20°  latitude  by  20°  longitude  area  and  was  remapped 
into  a  mercator  projection  using  the  program  REMAP  (Fig.  3.2). 

The  remapped  image  contained  a  large  amount  of  cloud  imformation,  unfor- 
tunately too  much  for  the  DISPLAY  software.  To  reduce  the  number  of  clouds,  only 
those  with  a  pixel  value  of  160  or  greater  (  -53°  C  )  were  selected.  The  hook  shape  of 
the  occluded  low  in  Figs.  3.1  and  3.2  is  still  clear  after  the  majority  of  clouds  are 
removed  by  setting  the  pixel  threshold.  The  resulting  image  and  500  mb  contours, 
overlayed  on  the  remapped  satellite  image,  is  presented  in  Fig.  3.3. 

2.  Atmospheric  Parameter  Data 
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Fig.  3.1     IR  GOES  Satellite  Image  (Raw) 
2300  GMT  15  May  1986. 
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Fig.  3.2     Mercator  Projected  Satellite  Image 
2300  GMT  15  May  1986. 
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Fig.  3.3     Mercator  Remapped  Satellite  Image, 
500  mb  Contours,  and  Enhanced  Pixels  (  Pixel  Value  ^  160). 
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The  two  dimensional  gridded  parameter  file  contours  are  presented  in  relation 
to  the  satellite  image  and  geopolitical  boundaries  in  Figs.   3.4,  3.5,  3.6,  and  3.7. 

The  height  range  of  the  500  mb  pressure  surface  (Fig.  3.4)  is  from  5640  to 
5880  meters  (m).  The  highest  point  on  that  surface  is  in  the  southwest  corner  and  the 
lowest  point  is  in  the  northwest  corner.  This  height  range,  240  m,  is  quite  small 
compared  to  the  surface  area  which  is  measured  in  hundreds  of  kilometers. 

The  320°  K  and  330°  K  constant  potential  temperature  surfaces  (Figs.  3.5 
through  3.7)  show  a  much  larger  vertical  variation  than  the  500  mb  pressure  surface. 
The  320°  K  surface  ranges  from  a  high  of  9498  m  in  the  northwest  corner,  to  a  low  of 
3165  m  in  the  southwest  corner.  The  330°  K  surface  ranges  from  a  high  of  11,150  m  in 
the  northwest  corner  and  to  a  low  of  6589  m  in  the  southeast.  Finally,  Fig.  3.7  is 
presented,  showing  both  the  320°  K  and  330°  K  surfaces.  Using  potential  temperature 
surfaces,  the  relationship  between  frontal  boundaries  and  cloudiness  can  be  examined. 

The  parameter  height  fields  were  then  processed  by  the  programs  GRIDIMG 
and  REMAP  to  produce  remapped  continuous  tone  images  for  use  as  inputs  to  the 
MOVIE2  programs  DEFSCT,  INTCLD,  CLDGEN  and  MERGER.  The  result  were 
files  that  described  the  geometry  of  the  3-D  structures  in  a  form  usable  as  inputs  to  the 
DISPLAY  program. 

E.       DISPLAY  PROCESSING 

To  demonstrate  the  information  conveyed  by  the  3-D  atmospheric  displays,  a 
series  of  frames  were  produced  at  various  angles  and  apparent  distances.  The  vertical 
height  scale  in  each  3-D  display  was  exaggerated  by  a  factor  of  50,  to  emphasize  the 
vertical  variations.  Each  figure  produced  consists  of  a  continuous  tone  image  and 
overlaying  line  graphic  display,  both  from  the  same  perspective.  The  overlaying 
process  helps  convey  additional  information  on  the  spatial  variation  of  the  structure 
being  viewed. 

In  each  image,  the  location  of  the  light  source  is  selected  to  be  either  at  infinity, 
in  a  given  direction,  or  at  the  eye  of  the  observer.  This  is  entered  during  execution  of 
the  DISPLAY  program  using  the  command  LIGHT.  The  user  describes  the  color  of 
each  part  withing  the  frame  using  the  COLOR  command  (Christiansen  and 
Stephenson,  1981b).  The  intensity  of  the  reflected  light  from  a  panel,  on  a  selected 
part,  is  dependent  on  the  orientation  of  the  panel.  With  a  variation  in  the  side  struc- 
ture of  clouds  or  a  height  change  on  a  data  surface,  the  result  is  light  returning  to  the 
viewer  at  different  intensities,  depending  on  the  angle  of  a  panel. 
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Fig.  3.4    500  mb  Pressure  Contours 
and  IR  Clouds  (  Pixels  2t  160  ). 
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Fig.  3.5     320°  K  Potential  Temperature  Contours 
and  IR  Clouds  (  Pixels  >  160  ). 
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Fig.  3.6    330°  K  Potential  Temperature  Contours 
and  IR  Clouds  (  Pixels  >  160  ). 
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Fig.  3.7     320°  K  and  330°  K  Contours 
and  IR  Clouds  (  Pixels  >  160  ). 
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Although  a  given  data  surface  may  be  defined  by  a  single  color,  varying  bright- 
nesses are  reflected  back.  The  result  of  these  shade  changes  is  that  information  about 
the  structure  of  the  surface  is  conveyed  to  the  viewer  which  might  otherwise  not  be 
apparent.  To  reenforce  the  changes  in  structure  and  depth,  line  graphic  images  were 
overlaid  on  the  continuous  tone  images.  This  gives  an  added  sense  of  depth  due  to  the 
divergent  and  receding  line  segments. 

In  producing  each  of  the  3-D  displays,  the  cloud  field  is  averaged.  Clouds  with 
pixel  values  of  160  or  greater  (Figs.  3.3  and  3.7)  are  smoothed.  The  result  is  a  loss  of 
the  smaller  cloud  tops  in  the  central  and  southern  portions  of  the  area.  However,  the 
hook  shape  is  still  evident  in  the  3-D  displays  (Figs.  3.10  and  3.13). 

Fig.  3.8  is  a  3-D  display  of  the  information  contained  in  Fig.  3.4.  In  this  display 
the  data  structures  (frame  parts)  are  viewed  from  the  southern  edge.  This  display 
places  in  proper  perspective  the  500  mb  pressure  surface,  clouds,  and  topography. 

To  examine  the  topographic  and  parameter  surface  variations  more  closely,  the 
apparent  distance  to  the  frame  was  reduced  by  one  half  (Fig.  3.9).  In  this  view,  the 
topographic  variations  are  slightly  more  distinct.  The  500  mb  surface  variations  are 
basically  indistinguishable,  even  with  the  factor  of  50  applied  to  the  vertical.  To  obtain 
a  view  which  conveys  more  information  about  the  500  mb  surface,  an  oblique  view  was 
produced,  at  the  original  distance  to  the  frame  (Fig.  3.10).  However,  here  again  the 
vertical  changes  in  the  500  mb  surface  can  barely  be  detected. 

As  described  earlier,  the  constant  potential  temperature  surfaces  have  consider- 
ably more  vertical  variation  than  the  pressure  surface.  Figures  3.11  through  3.13  show 
the  same  perspective  changes  as  with  the  500  mb  series.  A  comparison  of  Figs.  3.7  and 
3.13  illustrates  the  significant  increase  in  visual  clarity  and  understanding  of  the  rela- 
tionship between  the  two  potential  temperature  surfaces.  Fig.  3.7  must  be  examined 
carefully  and  for  some  time  before  their  relationship  can  be  established. 
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Fig.  3.8     3-D  500  mb  Surface,  Clouds  and  Topography 
View  From  the  Southern  Edge  of  the  Region. 
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Fig.  3.9    Zoomed  500  mb  Surface,  Clouds  and  Topography 
View  From  the  Southern  Edge  of  the  Region. 
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Fig.  3.10    Oblique  500  mb  Surface,  Clouds  and  Topography 
View  From  the  Southwest  Corner  of  the  Region. 
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Fig.  3.11     3-D  320°  K  and  330°  K  Surfaces 
View  From  the  Southern  Edge  of  the  Region. 
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Fig.  3.12    Zoomed  Potential  Temperature  Surfaces 
View  From  the  Southern  Edge  of  the  Region. 
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Fig.  3.13    Oblique  Potential  Temperature  Surfaces 
View  From  the  Southwest  Corner  of  the  Region. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  series  of  3-D  displays  of  the  atmospheric  data,  used  in  conjunction  with 
conventional  2-D  displays,  provide  another  tool  for  use  ■  by  meteorologists  in  their 
attempt  to  understand  weather  systems.  These  have  application  for  the  display  of 
pressure  surfaces.  Those  displays,  with  a  series  of  vertical  height  exaggerations,  could 
be  instructive  to  students  learning  about  the  structure  of  the  atmosphere.  The  infor- 
mation conveyed  in  an  operational  setting  would  be  useful  as  an  aid  to  flight  briefings. 
Aircrews  planning  flights  at  constant  pressure  levels  would  see  the  relationships 
between  their  flight  path  and  cloud  patterns  at  that  flight  level. 

Some  data  surfaces  vary  greatly  in  height  over  relatively  small  horizontal 
distances.  The  display  of  potential  temperature  surfaces  shows  the  possible  application 
of  such  models  to  atmospheric  studies.  This  depends  on  the  type  of  surface  and  inten- 
sity of  the  weather  situation  in  which  it  occurs.  Another  type  of  surface  that  varies 
greatly  in  height  are  constant  equivalent  potential  temperature  surfaces.  These  include 
atmospheric  moisture  in  their  computation  and  are  useful  in  locating  frontal  bound- 
aries at  various  points  both  vertically  and  horizontally.  However,  at  the  present  time 
MOVIE2  programs  cannot  adequately  deal  with  points  on  the  surface  topography 
intersected  by  such  isentropic  surfaces. 

Numerical  weather  prediction  (NWP)  models,  whose  output  is  the  basis  of  much 
of  present-day  weather  forecasting,  predict  a  variety  of  surfaces  which  could  be  viewed 
from  the  three-dimensional  perspective.  This  would  convey  added  information  about 
the  behavior  and  structure  of  the  various  levels  within  the  models.  Both  research  and 
operational  meteorologists  would  find  this  application  useful. 

B.  RECOMMENDATIONS 

After  having  developed  the  procedure  for  including  data  surfaces  in  these  3-D 
images  and  expanding  the  technique  to  a  synoptic  scale,  it  has  become  apparent  that 
many  improvements  are  possible.  The  current  sequence  and  time  for  construction  of 
these  images  needs  to  be  improved.  At  present,  a  number  of  steps  performed  with 
much  human  intervention  should  be  automated.  This  could  be  accomplished  by  estab- 
lishing default  or  standardized  processing  parameters. 
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Prior  to  actual  use  in  the  DISPLAY  software,  the  processing  could  occur  as  a 
result  of  program  executions  from  previously  constructed  files.  These  files  would 
control  program  execution,  providing  the  information  on  input/output  file  names  and 
settings  required  by  each  program  to  produce  the  file  for  input  to  DISPLAY.  This 
would  make  it  easier  for  users  to  employ  the  3-D  technique  in  case  and  model  studies. 

Within  the  program  DISPLAY,  several  modifications  are  recommended.  First, 
the  axis  triad,  showing  the  location  of  the  x,  y  and  z  axes,  should  include  a  vertical 
height  scale  along  the  z-axis.  This  would  greatly  improve  the  user's  sense  of  height  to 
width  proportions.  Second,  an  overlay  of  geopolitical  boundaries  within  the  3-D  image 
would  aid  in  properly  placing  features  with  respect  to  surface  locations.  Third,  an 
option  could  be  included  in  the  DISPLAY  software  allowing  simultaneous  construction 
and  display  of  the  line  graphic  and  continuous  tone  images.  As  can  be  seen  in  the 
figures  above,  the  combination  of  line  and  continuous  tone  images  provides  more 
information  on  the  structure  of  clouds  and  surfaces  than  either  alone.  Fourth,  an 
option  to  send  constructed  images  to  disk  or  tape  files  from  within  DISPLAY  should 
be  developed.  At  the  current  time,  the  user  must  exit  the  DISPLAY  program  abruptly 
to  ensure  that  the  3-D  images  are  not  lost,  as  when  normal  exit  from  DISPLAY  is 
used.  Saving  the  3-D  images  allows  very  rapid  retrieval  and  would  be  easier  to  use  in  a 
looping  process  for  time  sequence  examinations. 

The  CSU  generated  software  allows  for  only  three  layers  of  data.  That  number 
could  be  raised,  as  long  as  the  total  number  of  parts  in  the  DISPLAY  frame  do  not 
exceed  27.  A  capability  to  process  atmospheric  data  surfaces  as  unique  parts,  as 
opposed  to  modified  cloud  structures  should  be  added  to  the  programs  DEFSCT, 
INTCLD,  and  CLDGEN.  The  program  CLDGEN  should  be  modified  to  set  the  base 
of  pressure  and  other  surfaces  to  equal  the  topography  gray  shade  at  intersection 
points.  This  would  allow  the  display  of  any  atmospheric  surface,  especially  those 
whose  computed  elevations  may  be  below  the  surface.  The  actual  height  would  still  be 
available  in  the  grided  data  file  and  the  displayed  3-D  image  would  show  intersection 
locations. 

The  hardware  that  presently  produces  the  3-D  frames,  a  VAX  11/780  computer, 
is  a  time  shared  device.  Each  user  has  a  "virtual  machine",  seeing  what  appears  to  be  a 
computer  devoted  solely  to  that  individuals  task.  In  fact,  it  is  working  on  a  variety  of 
tasks  for  different  users,  allocating  its  attention  from  one  to  another  task.  The 
construction  of  3-D  frames  requires  the  use  of  a  large  amount  of  the  VAX  1 1/780 
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memory.  With  the  present  system  configuration  it  can  take  as  little  as  two  or  as  long 
as  long  as  thirty  minutes  to  construct  a  single  3-D  frame.  The  recommended  solution 
is  to  move  the  DISPLAY  software,  at  least,  onto  an  array  processor.  That  device 
could  do  the  complex  matrix  operations  involved  in  image  manipulations  in  much  less 
time. 

An  increase  in  the  resolution  of  the  3-D  displays  will  allow  active  use  by  NWP 
modelers,  researchers  studying  significant  weather  events  and  meteorology  students. 
With  some  modification,  oceanographic  data  fields  could  be  used  as  input  to  the  3-D 
display  technique. 

Sophisticated  3-D  display  hardware  is  becomming  available  for  researchers. 
Using  the  3-D  display  data  files  as  input,  devices  will  be  capable  of  letting  users  take 
real-time  animated  tours  through  the  displays.  This  should  lead  to  new  techniques  for 
the  analysis  of  data  that  previously  has  been  represented  in  only  two  dimensions. 
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APPENDIX  A 
SAMPLE  PROCESSING  STEPS 

For  those  users  who  would  like  to  produce  3-D  images,  the  following  lists  of 
suggested  processing  steps  is  offered  to  assist.  It  is  assumed  that  the  reader  is  familiar 
with  operating  on  the  VAX  1 1/780  computer,  using  the  COMTAL  display,  and  using 
GEMPAK. 

1.        CLOUDS  AND  TOPOGRAPHY 

•  Extract  the  topography  tape  records  using  the  MOVIE2  program  DMADAT. 
The  program  DMA512,  also  in  the  MOV1E2  package,  does  not  currently  write 
out  the  correct  header  record  and  should  not  be  used.  The  output  from 
DMADAT  is  a  file  of  your  name  choice  with  the  extension  ".STI". 

•  Construct  the  topography  model  using  the  output  file  from  step  1  as  input  to 
TOPGEN.  Note:  Any  vertical  height  exaggeration  factor  applied  here  will  be 
also  applied  to  subsequently  developed  cloudor  atmospheric  data  plane.  Output 
is  a  file  with  extension  ".STM '. 

•  Run  the  program  C '  r^CT  using  the  ".STI"  topography  file  and  a  cloud  image 
file  ".IMG"  as  input.  Both  a  visible  and  an  IR  cloud  image  file  can  be  used,  or 
an  IR  image  alone. 

•  In  keeping  track  of  the  many  files  produced  in  this  overall  sequence,  it  is  useful  to 
use  the  same  filename  for  each  file.  The  distinction  between  various  file  types  is 
made  in  the  extension  which  the  MOVIE2  software  assigns  to  the  output  files. 
Also,  the  topography  ".STI"  filename  and  the  ".CTM"  cloud  file  name  must 
match  for  them  to  oe  merged  together  later. 

•  The  side  length,  smoothing  factor  and  x,v  coordinates  of  the  top  left  corner  of  the 
TOPGEN  and  bEFSCT  files  must  match. 

•  Run  INTCLD  using  the  previous  steps'  output  files.  Make  sure  that  the 
cloud/no  cloud  threshold  input  value  is  a  gray  shade  appropriate  for  the  clouds 
you  wish  to  display.  If  the  threshold  is  set  too  low,  this  may  cause  too  many 
parts  to  be  included  in  the  3-D  model  (27  is  the  maximum  number).  Set  the 
cloud  bases  to  a  level  determined  by  examination  of  a  representative  rawinsonde 
sounding  for  the  area. 

•  Run  CLDGEN.   No  prompts  should  appear  in  this  step. 

•  Run  MERGER.  No  prompts  at  this  step  either.  Note  the  number  of  clouds 
merged.  This  information  is  used  later  in  DISPLAY  to  identify  parts  of  the  3-D 
frame  on  which  to  operate.    ' 

•  Run  the  DISPLAY  program  using  the  command  "R  MOVIE2". 

•  Once  in  DISPLAY,  the  user  is  prompted  for  an  input  file  name.  The  response 
can  be  either  the  topography  ".STM'  file  to  a  cloud  topography  model  ".CTM" 
file. 

•  Within  the  DISPLAY  software,  commands  can  be  referenced  in  the 
MOVIE. BYU  Training  Manual.  However,  examples  of  the  responses  to  these 
commands  are  difficult  to  locate.  When  producing  continuous  tone  color  frames, 
the  user  must  enter  the  COLOR  and  LIGHT  commands  and  respond  to  a  series 
of  prompts.  An  example  of  the  response  to  the  COLOR  command  prompts 
appears  below.    User  responses  do  not  include  the  quote  marks. 
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—Background— (Red, Green, Blue)—  >  >  >"0,.7,.7" 
— Parts--- >  >  >\l  .6A.6" 

>  >  >  "2,3  0,.9,(3" 

>  >  >  "6,6  0,0,.9" 

>  >  >"1,10,.1,.1" 

•  After  the  next  >  >  >  prompt,  the  user  enters  a  null  entry. 

— Standard  Fringe— >  >  "Yes" 
—Symmetric >  >  "Yes" 

•  Typical  responses  to  the  LIGHT  command  might  be  as  follows. 

—  Light  Source  Number - >  >  >"\" 

—Light  Source  at  the  eye  of  the  observer? — >  >  >  "Yes" 

—Light  Source  Intensity >  >  >".8" 

—Normal  Parts  Exponent >  >  >"l,7  .6" 

—Highlight  for  parts. ..and  exponent >  >  >  "1,7  .6,6" 

•  After  the  last  >  >  >  the  user  sends  a  null  entry. 

CLOUDS,  TOPOGRAPHY,  DATA  PLANES 

•  Produce  gridded  data  sets  for  this  procedure  in  GEMPAK.  The  area  enclosed  by 
the  boundaries  of  the  gridded  data  set  should  be  about  5  degreees  larger  on  each 
side  than  the  satellite  and  topography  images.  This  allow  for  the  remapping 
process. 

•  Make  note  of  the  maximum  and  minimum  heights  in  the  grided  data  file.  Match 
these  heights  with  the  corresponding  temperature  on  a  representative  rawinsonde 
sounding.   Convert  these  temperatures  into  IR  pixel  values. 

•  Create  navigation  files  for  the  data  image  using  the  program  NAVGEN. 

•  Use  the  data  plane  image  and  the  navigation  files  as  input  to  the  GEMPAK 
program  REMAP.  The  latitude  and  longitude  of  the  remapped  image  bounds 
should  match  those  in  the  topography  andcloud  image  files. 

•  Run  the  output  mercator  projected  image  through  DEFSCT,  INTCLD,  and 
CLDGEN  as  previouslv  described.  In  INTCLD  enter  a  base  height  of  "-999". 
Delete  the  DEFSCT.DAT  file  ater  each  program  is  executed. 

Rename  the  .BLO  and  .LOW  files  output  from  CLDGEN  as  the  ".BMI"  and 
".MID",  or  ".BHI"  and  ".HGH"  file,  respectively. 

Run  MERGER. 

Run  DISPLAY  as  previously  described. 


• 
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APPENDIX  B 
PREPROCESSING  OF  SATELLITE  DATA 

The  satellite  data  originating  at  NEPRF  is  on  800  bpi- tapes.    In  order  to  read  the 

tape  files  onto  the  VAX  computer,  they  must  be  converted  into  1600  bpi.   This  process 

occurs  on  the  NPS   UNIX  computer.    The  list  below  describes  steps  used  in  this 

process  that  may  be  helpful. 

Log  onto  the  UNIX  computer  using  the  "GUEST"  account.    The  password  may 
be  obtained  from  the  UNIX  system  manager. 

•  Create  a  directory  to  temporarilv  contain  the  satellite  and  navigation  files  using 
the  following  command:  MKD1R  TEMP".  Any  directory  name  can  be  used  in 
place  of  the   TEMP"  entry. 

•  Change  to  the  "TEMP"  directory  by  issuing  the  command  "CD  TEMP". 

•  After  mounting  the  tape  on  the  UNIX  tape  drive  enter  the  following  command: 
"dd  if=/dev/rmt4  of=  fname".  "fname"  should  be  replaced  with  a  unique  name 
for  each  file  read  off  the  tape.  The  above  command  must  be  issued  for  each  tape 
file.  The  UNIX  system  will  prompt  the  user  for  a  new  filename  after  each 
previous  file  has  been  read  in.  This  process  reads  in  the  800  bpi  image  and  navi- 
gation files  to  the  files  named  by  the  user,  in  the  directory  created  at  the  begin- 
ning. Note:  UNIX  interprets  upper  and  lower  case  letters  differently.  User's 
should  enter  the  commands  in  lower  case. 

•  Rewind  the  tape. 

•  Make  sure  the  tape  "write  ring"  is  on  the  data  tape. 

•  Read  each  disk  file  back  onto  the  tape  using  the  command  "dd  if=  fname 
of=/dev/rmtl2".  This  writes  the  files  back  to  tape  in  1600  bpi  form.  Similar  to 
before,  '  fname"  should  be  substituted  with  the  name  of  the  image  or  navigation 
file  to  be  written  out  to  tape. 

•  Rewind. the  tape  and  log  off  the  UNIX  computer  using  the  command  "logout". 

•  Log  onto  the  VAX  computer. 

•  Enter  the  command  "ALLOCATE  MFA0:"  to  reserve  the  VAX  tape  drive  for 
use.  This  ensures  that  another  user  will  not  accidentally  write  to  the  satellite 
data  tape  and  erase  files. 

•  Assign  the  MFA0  device  as  a  tape  reading  device  using  the  command  "ASSIGN 
MFA0:  TAPE". 

•  Enter  "MOUNT/FOR/DEN  =  1600  MFA0:". 

•  Enter  "RUN  TAPES:  MAGTODSK".  This  command  executes  a  VAX  utility  to 
read  magnetic  tape  disk  files  sequentially  to  disk  files. 

•  In  response  to  the  record  length  prompt,  the  user  should  enter  a  "N". 

•  The  input  device  name  is  "MF0:".  Note  that  there  is  no  "A"  as  in  "MFA0". 

•  Specify  the  desired  output  disk  file  names  when  prompted.  This  must  be  done  for 
each  tape  file. 
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•  Enter  the  command  "DISMOUNT  MFAO:"  and  then  "DEALLOCATE  MFAO:". 
This  first  rewinds  the  tape  and  then  makes  the  tape  drive  available  for  other 
users. 

•  After  the  images  and  navigation  files  are  determined  to  be  correct.  The  user  must 
log  back  onto  the  UNIX  computer  and  delete  the  directory  containing  the 
temporary  image  and  navigation  files. 
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